A1 = readmatrix('NiI2_CPGE_680nm_-24V_30K_2.txt')
A2 = readmatrix('NiI2_CPGE_680nm_-24V_40K.txt')
A3 = readmatrix('NiI2_CPGE_680nm_-24V_50K.txt')
A4 = readmatrix('NiI2_CPGE_680nm_-24V_60K.txt')

x1 = A1(:,1)
y1 = A1(:,2)*1e12

x2 = A2(:,1)
y2 = A2(:,2)*1e12-0.1

x3 = A3(:,1)
y3 = A3(:,2)*1e12-0.2

x4 = A4(21:60,1)-10
y4 = A4(21:60,2)*1e12

hold on
xlim([11, 50])
plot(x1,y1,'Linewidth',1,'Color',[0.2 0 0.8])
plot(x2,y2,'Linewidth',1,'Color',[0.4 0 0.6])
plot(x3,y3,'Linewidth',1,'Color',[0.6 0 0.4])
plot(x4,y4,'Linewidth',1,'Color',[0.8 0 0.2])
xlabel('Polarization')
ylabel('Photocurrent [pA]')
ylim([0,2.1])
yticks(0:1:2)
%set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
hold off

set(gcf,'Position',[0, 0, 190, 190])